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Abstract —This paper addresses the problem of deriving the 
asymptotic performance of adaptive Low Rank (LR) filters used 
in target detection embedded in a disturbance composed of a 
LR Gaussian noise plus a white Gaussian noise. In this context, 
we use the Signal to Interference to Noise Ratio (SINR) loss 
as performance measure which is a function of the estimated 
projector onto the LR noise subspace. However, although the 
SINR loss can be determined through Monte-Carlo simulations 
or real data, this process remains quite time consuming. Thus, 
this paper proposes to predict the SINR loss behavior in order to 
not depend on the data anymore and be quicker. To derive this 
theoretical result, previous works used a restrictive hypothesis 
assuming that the target is orthogonal to the LR noise. In 
this paper, we propose to derive this theoretical performance 
by relaxing this hypothesis and using Random Matrix Theory 
(RMT) tools. These tools will be used to present the convergences 
of simple quadratic forms and perform new RMT convergences of 
structured quadratic forms and SINR loss in the large dimensional 
regime, i.e. the size and the number of the data tend to infinity at 
the same rate. We show through simulations the interest of our 
approach compared to the previous works when the restrictive 
hypothesis is no longer verified. 

Index Terms —Low Rank SINR loss, Random Matrix Theory, 
Adaptive Filtering, Quadratic Forms convergence, Spiked model 


I. Introduction 

In array processing, the covariance matrix R of the data 
is widely involved for main applications as filtering [1], [2], 
radar/sonar detection [3] or localization [4], [5]. However, 
when the disturbance in the data is composed of the sum 
of a Low Rank (LR) correlated noise and a White Gaussian 
Noise (WGN), the covariance matrix is often replaced by the 
projector onto the LR noise subspace II C [6]—[9]. In practice, 
the projector onto the LR noise subspace (and the covariance 
matrix) is generally unknown and an estimate is consequently 
required to perform the different processing. This estimation 
procedure is based on the so-called secondary data assumed 
to be independent and to share the same distribution. Then, 
the true projector is replaced by the estimated one in order to 
obtain an adaptive processing. An important issue is then to 
derive the theoretical performances of the adaptive processing 
as a function of the number of secondary data K. The 
processing based on the covariance matrix has been widely 
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studied and led to many theoretical results in filtering [1] 
and detection [10] — [13]. For example, for classical adaptive 
processing, K = 2m secondary data (where m is the data 
size) are required to ensure good performance of the adaptive 
filtering, i.e. a 3dB loss of the output Signal to Interference 
plus Noise Ratio (SINR) compared to optimal filtering [1], 
For LR processing, some results has been obtained especially 
in filtering [6], [14]—[16] and localization [17], Similarly, in 
LR filtering, the number K of secondary data required to 
ensure good performance of the adaptive filtering is equal to 
2r (where r is the rank of the LR noise subspace) [6], [14], 

These last results are obtained from the theoretical study 
of the Signal to SINR loss. More precisely, in [14], [16], the 
derivation of the theoretical results is based on the hypoth¬ 
esis that the steering vector is orthogonal to the LR noise 
subspace. Nevertheless, even if the result seems to be close 
to the simulated one when the hypothesis is no longer valid 
anymore [18], it is impossible with traditional techniques 
°f [14], [16] to obtain a theoretical performance as a function 
of the distance between the steering vector and the LR noise 
subspace. Since, in practice, this dependence is essential to 
predict the performance of the adaptive filtering, we propose in 
this paper to derive the theoretical SINR loss, for a disturbance 
composed of a LR noise and a WGN, as a function of K and 
the distance between the steering vector and the LR noise 
subspace. The proposed approach is based on the study of the 
SINR loss structure. 

The SINR loss (resp. LR SINR loss) is composed of 
a simple Quadratic Form (QF) in the numerator, .s^R Ls^ 
(resp. sf 1 IQ^) and a structured QF in the denominator 
sf r R _1 RR' 1 S 2 (resp. RrQs 2 ). These recent years, 

the simple QFs (numerator) have been broadly studied [19]— 
[22] using Random Matrix Theory (RMT) tools contrary to 
structured QFs (denominator). RMT tools have also been used 
in array processing to improve the MUSIC algorithm [23], [24] 
and in matched subspace detection [25], [26] where the rank r 
is unknown. The principle is to examine the spectral behavior 
of R by RMT to obtain their convergences, performances and 
asymptotic distribution when K tends to infinity and when 
both the data size to and K tend to infinity at the same 
ratio, i.e. m/K -A c g]0, +oo), for different models of R 
of the observed data as in [19], [20], [23], [22] and [21], 
Therefore, inspired by these works, we propose in this paper to 
study the convergences of the structured QFs s^R. _1 RR 
(resp. .s^rQRrQ^): when 1) K — > oo with a fixed to 
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and when 2) to, K —> oo at the same ratio under the most 
appropriated model for our data and with the rank assumed 
to be known. From [27], [28], the spiked model has proved 
to be the more appropriated one to our knowledge. This 
model, introduced by [29] (also studied in [30], [31] from 
an eigenvector point of view) considers that the multiplicity 
of the eigenvalues corresponding to the signal (the LR noise 
for us) is fixed for all to and leads to the SPIKE-MUSIC 
estimator [32] of .s^ILs^. Then, the new results are validated 
through numerical simulations. From these new theoretical 
convergences, the paper derives the convergence of the SINR 
loss for both adaptive filters (the classical and the LR one). 
The new theoretical SINR losses depend on the number of 
secondary data I\ but also on the distance between the steering 
vector and the LR noise subspace. This work is partially 
related to those of [33], [34] and [35] which uses the RMT 
tools to derive the theoretical SINR loss in a full rank context 
(previously defined as classical). 

Finally, these theoretical SINR losses are validated in a 
jamming application context where the purpose is to detect a 
target thanks to a Uniform Linear Antenna (ULA) composed 
of to sensors despite the presence of jamming. The response 
of the jamming is composed of signals similar to the target 
response. This problem is very similar to the well-known 
Space Time Adaptive Processing (STAP) introduced in [2], 
Results show the interest of our approach with respect to other 
theoretical results [6], [14]—[16] in particular when the target 
is close to the jammer. 

The paper is organized as follows. Section II presents the 
received data model, the adaptive filters and the corresponding 
SINR losses. Section III summarizes the existing studies 
on the simple QFs .s^R.s^ and .s^ll.sa, and exposes the 
covariance matrix model, the spiked model. Section IV gives 
the theoretical contribution the paper with the convergences 
of the structured QFs and and 

the convergences of the SINR losses. The results are finally 
applied on a jamming application in Section V. 

Notations: The following conventions are adopted. An 
italic letter stands for a scalar quantity, boldface lower¬ 
case (uppercase) characters stand for vectors (matrices) and 
{.) H stands for the conjugate transpose. Ijv is the N x 
N identity matrix, tr(.) denotes the trace operator and 
diag(.) denotes the diagonalization operator such as (A)^ = 
(diag(a)) M = (a)* and (A)^- = 0 if i ± j. # {A} denotes 
the cardinality of the set A. [a, 6 ] is the set defined by 
{cc £ Z : a ^ x < 6 , V(a, b) £ Z 2 }. <D nX N is a nx N matrix 
full of 0. The abbreviations iid and a.s. stem for independent 
and identically distributed and almost surely respectively. 

II. Problem statement 


is an unknown deterministic vector containing the different 
parameters of the target (e.g. the localization, the velocity, 
the Angle of Arrival (AoA), etc.). In the remainder of the 
article, in order to simplify the notations, © will be omitted 
of the steering vector which will simply be denoted as a. If 
necessary, the original notation will be taken. 

This section will first introduces the data model. Then, the 
filters, adaptive filters and the corresponding SINR loss, the 
quantity characterizing their performances, will be defined. 


A. Data model 

The observation vector can be written as: 


x = d + c + b 


(1) 


where c + b is the noise that has to be whitened, b £ 
C mxl ~ CAf(0,er 2 I m ) is an Additive WGN (AWGN) and 
c is a LR Gaussian noise c £ C mxl modeled by a zero- 
mean complex Gaussian vector with a normalized covariance 
matrix C (tr(C) = to), i.e. c ~ CAf(0, C). Consequently, 
the covariance matrix of the noise c + b can be written as 
R = C 4 - a 2 I m £ C mxm . Moreover, considering a LR 
Gaussian noise, one has rank (C) = r <C to and hence, the 
eigendecomposition of C is: 

r 

c=y: ii u i u ? (2) 

i= 1 

where 7 ,; and it,, i £ [l;r] are respectively the non-zero 
eigenvalues and the associated eigenvectors of C, unknown 
in practice. This leads to: 

m 

R = XiUiuf (3) 

i=l 


where A* and Ui, i £ [l,m] are respectively the eigenvalues 
and the associated eigenvectors of R with Ai = 71 + er 2 > 
• • • > X r = 7 r + a 2 > A r+ i = ■■■ = X m = a 2 . Then, the 
projector onto the LR Gaussian noise subspace II C and the 
projector onto the orthogonal subspace to the LR Gaussian 
noise subspace II ( L are defined as follows: 


n c = £l=i Uiuf 

n f A = i m - n c = £™ r+1 u lU » 


(4) 


However, in practice, the covariance matrix R of the noise 
is unknown. Consequently, it is traditionally estimated with the 
Sample Covariance Matrix (SCM) which is computed from 
K iid secondary data Xk = c k + b k , k £ [1, A'], and can be 
written as: 


^ K m 

R = — x k x% = 

k—1 i—1 


(5) 


The aim of the problem is to filter the received observation 
vector x £ C mxl in order to whiten the noise without 
mitigating an eventual complex signal of interest d (typically 
a target in radar processing). In this paper, d will be a target 
response and is equal to aa(&) where a is an unknown 
complex deterministic parameter (generally corresponding to 
the target amplitude), a(0) is the steering vector and 0 


where Ai and u, t , i £ [l,m] are respectively the eigenvalues 
and the eigenvectors of R with Ai A A 9 A ■■■ A A m , c k ~ 
CN{0 ,C) and b k ~ CA"(0- <r 2 I m ). Finally, the traditional 
estimated projectors estimated from the SCM are: 


ri c = EL-i un if 

n c ^ = i m -n c = £™ r+ i^f, 


(6) 
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B. Adaptive filters 

A filtering preprocessing on the observation vector x 
(Eq.(l)) is first done with the filter w in order to whiten the 
received signal p = w H x. The filter maximizing the SINR is 
given by: 

w opt = R : a (7) 

Since, in practice, the covariance matrix R of the noise is 
unknown, the estimated optimal filter or adaptive filter (sub- 
optimal) is: 

w = R _1 a (8) 

In the case where one would benefit of the LR structure 
of the noise, one should use the optimal LR filter, based on 
the fact that II t 4 is the best rank r approximation of R 1 , 
which is defined by [6]: 

™lr = n^a (9) 

As, in practice, the projector is not known and is estimated 
from the SCM, the estimated optimal filter or adaptive filter 
(sub-optimal) is: 

i&lr = (10) 


C. SINR Loss 

Then, we define the SINR Loss. In order to characterize 
the performance of the estimated filters, the SINR loss com¬ 
pares the SINR at the output of the filter to the maximum 
SINR: 

SIN R ou t \w H d\ 2 

P SIN R max (w H Rw)(d H Li^ 1 d) 

\„ H d - 1„|2 

\a iv a 

(a ff R-'RR- 1 a)(a fl R- 1 a) 


(ID 

( 12 ) 


If w = w opt , the SINR loss is maximum and is equal to 1. 
When we consider the LR structure of the noise, the theoretical 
SINR loss can be written as: 


Plr = 


Knd\ 2 


(iu^Riulr) {d H R~ 1 d) 


a M ntra\ 


(a H n c i Rn c i a)(a Jf R- 1 a) 


(13) 


(14) 


Linally, the SINR loss corresponding to the adaptive filter in 
Eq.(10) is defined from Eq.(14) as: 


Plr — PLR.In-=nj- (15) 

Since we are interested in the performance of the filters, 
we would like to obtain the theoretical behavior of the SINR 
losses. Some asymptotic studies on the SINR loss in LR 
Gaussian context have already been done [14], [16]. In [14], 
[16], the theoretical result is derived by using the assumption 
that the LR noise is orthogonal to the steering vector and, in 
this case, [16] obtained an approximation of the expectation 
of the SINR loss /3lr- However, this assumption is not always 
verified, not very relevant and is a restrictive hypothesis in 
real cases. We consequently propose to relax it and study 
the convergence of the SINR loss using RMT tools through 
the study of the nominators and denominators. Indeed, one 
can already note that the numerators are simple QLs whose 


convergences were widely considered in RMT. However, the 
denominators contain more elaborated QLs which were not 
tackled in RMT yet and will be the object of Sec. IV. 

III. Random matrix theory tools 

This section is dedicated to the introduction of classical 
results from the RMT for the study of the convergence of QLs. 
This theory and the convergences are based on the behavior of 
the eigenvalues of the SCM when m, K —>■ oo at the same rate, 
i.e. m/K —> c £ ]0, + 00 ). In order to simplify the notations, 
we will abusively note c = m/K. 

The useful tools for the study of the eigenvalues behavior 
and the assumptions to the different convergences will be first 
presented. Secondly, the section will expose the data model, 
the spiked model [22]. Linally, the useful convergences of 
simple QLs (s A R _1 S 2 , •s( , II.S 2 ) will be introduced. 


A. Preliminaries 


The asymptotic behavior of the eigenvalues when 
to, K — > 00 at the same rate is described through the conver¬ 
gence of their associated empirical Cumulative Distribution 
Lunction (CDL) F m (x ) or their empirical Probability Density 
Lunction (PDL) f m (x ) 1 . The asymptotic PDL f m (x) will 
allow us to characterize the studied data model. The empirical 
CDL of the sample eigenvalues of R can be defined as: 

F m {x) = — #\k : A fc < z) (16) 

to 1 J 


However, in practice, the asymptotic characterization of F m ( x ) 
is too hard. Consequently, one prefers to study the convergence 
of the Stieltjes transform (SF [■]) of F m (x): 


bm (a) 


ST 

1 


E 

i= 1 


F m (x) 

1 


x — z 


-dF m (x) 


A i — z 


1 

= —tr 


(R - zl m ) 


-1 


(17) 

(18) 


with z £ C+ = {z £ C : 3 [2] >0} and which almost surely 
converges to b m (z). It is interesting to note that the PDL can 
thus be retrieve from the Stieltjes transform of its CDL: 


fm(x) 


lim —3 

S[z]->0 7T 



(19) 


with x £ R. In an other manner, the characterization of / m (x) 
(resp. fm(x )) can be obtained from b m (z ) (resp. b m (z)). Then, 
to prove the convergences, we assume the following standard 
hypotheses. 

(Asl) R has uniformly bounded spectral norm Vm £ N*, i.e. 

V* £ [1, to], A i < 00 . 

(As2) The vectors si, S 2 £ C mxl used in the QLs (here 
a(0) and x) have uniformly bounded Euclidean norm 

Vto £ N*. 

(As3) Let Y £ C mxA having iid entries with zero 
mean and unit variance, absolutely continuous and with 

Eflr/ijj 8 ] < 00 . 


^ne can show that under (Asl,As3-As5) described later, fm(x) a.s. 
converges towards a nonrandom PDF f(x) with a compact support. 
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(As4) Let Y £ €7" x h defined as in (As3), then its distribu¬ 
tion is invariant by left multiplication by a deterministic 
unitary matrix. Moreover, the eigenvalues empirical PDF 
of ^YY h a.s. converges to the Marcenko-Pastur distri¬ 
bution [36] with support [(1 — y/c) 2 , (1 + y/c) 2 ]. 

(As5) The maximum (resp. minimum) eigenvalue of f YY n 
a.s. tends to (1 + y/c) 2 (resp. to (1 — y/c) 2 ). 


B. Covariance matrix models and convergence of eigenvalues 

We first expose the considered data model and, then, the 
eigenvalues behavior of the SCM. The SCM can be written as 
R = with: 

X = R 1/2 Y = (I m + P) 1/2 Y (20) 

with X = [*!,••• ,xk\- The iid entries of Y follow the 
CAf( 0,1) distribution according to our data model in Sec. II. 
The complex normal distribution being a particular case of 
such distributions defined in (As3), the Y entries consequently 
verify it. Thus, the forthcoming convergences hold in the more 
general case defined by (As3). R 1 / 2 is the m x to Hermitian 
positive definite square root of the true covariance matrix. 
The matrix P is the rank r perturbation matrix and can be 
eigendecomposed as P = Uf2U H = YliLi with: 


WiIaCi 





( 21 ) 


with U = [Ui---Ujvj] and M the number of distinct 
eigenvalues of R. Moreover, Uj £ C mxlCi where /Cj is the 
multiplicity of Wj. Hence, the covariance matrix (3) can be 
rewritten as: 


M 

R = A * U * U f ( 22 ) 

2=1 

where Aof multiplicity AC,., and U, are the eigenvalues and 
the associated subspaces (concatenation of the ACj eigenvectors 
associated to Aj) of R respectively, with Ai = 1 + wi > • ■ ■ > 
Xm = 1 + Wjvr > 0 and = tn. The properties of the 

spiked model are the following: 

• 3n £ [1, M] such that co n = 0. 

• The multiplicity ACj is fixed Vz £ [1, MJ\n and does not 
increase with to, i.e. Ki/m —> 0 + , Vz£ [l,M]\n. 

m,K—foo 

Consequently, we have rank(fi) = X^:e|i Mj\n = r an d 
IC n = to — r. In other words, the model specifies that only 
a few eigenvalues are non-unit (and do not contribute to the 
noise unit-eigenvalues) and fixed. Consequently, X n = 1 is 
the eigenvalue of R corresponding to the white noise and the 
others correspond to the rank r perturbation. 

In our case (see Sec.II), we recall that the covariance 
matrix R can be written as in Eq.(3) and Eq.(22). More 
specifically, the noise component b corresponds to the white 
noise and its eigenvalue is A^j = 1 as, for simplicity pur¬ 
poses, we set < 7 " = 1. The r eigenvalues of the LR noise 
component c are strictly higher than 1. Thus, M = r + 1, 
Ai = 1 + W 1 > ■ • • > A M-l = 1 + W M -1 > A M = 1> ACj = 1 


is the multiplicity of Aj, Vz £ [1, r], and AC a7 = to. — r is the 
multiplicity of A a7 . 


M — l r 

R=A a7 U a 7 U^ + ^ AiUiUf=U r+ iU^ 1 +^Aiti i uf (23) 
2 = 1 2 = 1 

This model leads to a specific asymptotic eigenvalues PDF 
of R as detailed hereafter. The convergence of the eigenvalues 
is addressed through the convergence of the Stieltjes transform 
of the eigenvalues CDF. The asymptotic eigenvalue behavior of 
R for the spiked model was introduced by Johnstone [29] and 
its eigenvalue behavior was studied in [37]. In order to derive 
it, [37] exploited the specific expression given in Eq.(20). 
Then, [22] introduced the final assumption ( separation con¬ 
dition) under which the following convergences are given. 

(As 6 .S) The eigenvalues of P satisfy the separation condition, 
i.e. \uJi\ > yfc for all i £ [l,M]\n (z £ [1,7-J in our 
case). 

Thus, under (Asl-As5, As 6 .S), we have: 

fmix) —> f{x) (24) 


where f{x) is the Marcenko-Pastur law: 


/ 0 ) = < 


1 v^r 


2ttcx 

10 , 


if x = 0 and c > 1 
x)(x - A_|_), 

if x g]A_, A+[ 
otherwise 


(25) 


with A_ = (1 — y/c) 2 and A + = (1 + y/c) 2 . However, it is 
essential to note that, for all z £ [l,M]\n: 

\ a - s -, i , , 1 + W 2 

Aj'GAti -^ T i — 1 + Wj + C (26) 

m,K—> oo UJi 

where A ii is the set indexes corresponding to the j-th 
eigenvalue of R (for example A4 r +i = {r + 1, ■ • • , to} for 
A r +i). Two representations of f rn (x) for two different c and a 
sufficient large to are shown on Fig. 1 when the eigenvalues 
of R are 1, 2, 3, and 7 with the same multiplicity, where the 
eigenvalue 1 is the noise eigenvalue. One can observe that 
say (As 6 .S) is verified is equivalent to say that t„_. | > A+ 
and r„+i < A . In other words, all the sample eigenvalues 
corresponding to the non-unit eigenvalues of R, converge to a 
value Ti which is outside the support of the Marcenko-Pastur 
law (“asymptotic” PDF of the “unit” sample eigenvalues). As 
an illustration, one can notice that, in Fig. 1, for f m {x) plotted 
for c = 0.1, the separation condition is verified (uj i = 6, 
W 2 = 2 and ui?j = 1 are greater that y/c = 0.316) and the three 
non-unit eigenvalues are represented on the PDF and outside 
the support of the Marcenko-Pastur law by their respective 
limits ri = 7.116, T 2 = 3.15 and 73 = 2.2. On the contrary, for 
fm(x) plotted for c = 1.5, only the two greatest eigenvalues 
are represented on the PDF by their respective limits n = 8.75 
and T ‘2 = 5.25 while the separation condition is not verified 
for the eigenvalue A 3 = 2 (023 = 1 < yfc = 1.223). In this 
case, the sample eigenvalues corresponding to the eigenvalue 
A 3 = 2 belongs to the Marcenko-Pastur law. 
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Eigenvalues 

Fig. 1. PDF of the eigenvalues of the SCM with the spiked model when the 
eigenvalues of R. are 1, 2, 3, and 7 with the same multiplicity, where 1 is the 
noise eigenvalue. 


C. Convergence of simple quadratic forms 

Here, we compare the convergence of two QFs in two 
convergence regimes: when K — ¥ oo with a fixed m and when 
to, K —» oo at the same rate. 

We first present the useful convergences of simple QFs 
function of R. It is well known that, due to the strong law 
of large numbers, when K —t oo with a fixed m, R — > R 
a.s. [38]. Thus, 

sfR _1 S2 sfR -1 s 2 (27) 

K—foo 

m< oo 

Moreover, when to, K —>■ oo at the same rate [19], [39]: 

sfR _1 S 2 (l-c) -1 sfR -1 s 2 (28) 

m,K—f oo 

m / K —*c<oo 

The useful convergences of simple QFs function of IT, 1 are 
then presented. As R —> R a.s. when K —> oo with a fixed to, 
n ( x —> II, a.s. [19] in the same convergence regime. Thus: 

s fn x S2 sfn x s 2 (29) 

K^/- oo 


For the convergences in the large dimensional regime 
(to, K —> oo at the same rate), the convergences are presented 
under (Asl-As5) and the separation condition As6.S. [22] 
showed that, Vi € [1, M — 1]: 




«2 


1 — cue 




H t 


m,K—> oo 

7. / K ^ c < oo 


OjJa 


(30) 


with u>i = Xi — 1. A i is the i-th distinct eigenvalue of R. Let 

1 — CL >~ 2 

Xi = -' . . Thus, using the following relationship, 

1 + cw, 


n x = i ; 


M—l 

E UiU? = i 

i= 1 


E 

i—1 


one can deduce that with the spiked model and in the large 
dimensional regime: 

«fn x s 2 «fn x sS2 (32) 

m,K—f oo 

m / K —tc<oo 


with II X S = Yn=t 4’i u i u ? and 


4>i 


1, if i > r 

1 - Xi, if i sS r 


(33) 


Consequently, s^R _1 s 2 is consistent in the two con¬ 
vergence regimes and, although .s jQl, 1 .s 2 is consistent when 
K —> oc with a fixed to, it is no more consistent under the 
regime of interest i.e. when both to, K —> oc at the same rate. 


IV. New convergence results 


A. Convergence of structured quadratic forms 

In this section, the convergence of the structured QF 
function of fl, 1 is analyzed and results to Proposition 1. 

Proposition 1: Let B be a to x to deterministic complex 
matrix with a uniformly bounded spectral norm for all rn. 
Then, under (Asl-As5, As6.S) and the spiked model. 


sfn x Bn x s2 


sf nX s Bn x sS2 


(34) 


where II X S = 4’i u i u i I with defined by Eq.(33). 


Proof: See Appendix. 


Moreover, one can remark that if B = R, where R is 
the covariance matrix as defined in Eq.(5), the following 
convergence holds: 


sf n x Rn x s 2 


m,K—f oo 

n/K^c< oo 


sfn x s Rn c > 2 


(35) 


A visualization of the convergence of Eq.(34) in terms of Mean 
Squared Error (MSE) can be found in Fig. 2 when m, I\ —> oo 
at a fixed ratio. It is compared to the MSE corresponding to 
the following convergence when K —> oc with a fixed rn: 


(36) 


B. Convergence of SINR losses 

Now, we provide the convergences of the estimated SINR 
loss using the convergences previously presented and the 
following convergence. We recall that, as R -> R a.s. when 
K —> oo with a fixed to, one has: 

sfR _1 RR _1 s 2 sfR _1 s 2 (37) 

K—fOO 

m< oo 

Hence, when K —> oo with a fixed to and using Eq.(27), 
Eq.(37) and the continuous mapping theorem [38]: 

„ a.s. |q. h R~ 1 q| 2 

^ K-roc (a H H^ 1 a)(a H H~ 1 a) 


(31) 


(38) 
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Fig. 2. MSE over 10 3 iterations corresponding to Eq.(34) and Eq.(36) when 
the eigenvalues of Ft are 1, 21, 31, and 71 with the multiplicity m — 3, 1, 
1 and 1 respectively, c = 0.1, si = «2 are steering vectors of the LR noise 
component c and B = Ft. 


And, under (Asl-As5), when m, K —> oo at the same rate, 
from [33], we have: 


a.s.^ (1 - c)|a H R 1 a| 2 

m,K—f oo (a H R~ 1 a){a H 'R- 1 a) 

m/K —tc<oo 


(39) 


Thus, the estimated SINR loss p is consistent when K —> oo 
with m fixed and when m,K —>• oo at the same rate, up to 
an additive constant c. Consequently, RMT cannot help us to 
improve the estimation of the theoretical SINR loss. 

For the SINR loss corresponding to the adaptive LR filters, 
when K oo with a fixed m, using Eq.(29), Eq.(36) and the 
continuous mapping theorem theorem, we have: 

Plr -—> Plr (40) 

K —^oo 

m<oo 


where p\Ai is defined by Eq.(14). When m,K — > oo at the 
same rate, we obtain the following convergence: 

PLR Plr = Putln^n-L 7^ plr (41) 

m,K—> oo c c > b 

m / K —lc<oo 


of the jamming, c is composed of signals similar to the target 
response. In this section, except for the convergences when 
m, K —► oo at the same rate c, we choose m = 100 in 
order to have a large number for the data dimension. Even 
if, in some basic array processing applications, this number 
could seem significant, it actually became standard in many 
applications such as STAP [2], MIMO applications [40], [41], 
MIMO-STAP [40], etc. Here, 0 = 6 where 0 is the AoA. 
The jamming is composed of three synthetic targets with AoA 
—20°, 0° and 20° and wavelength l 0 = 0.667m. Thus, the 
jamming (LR noise) has a rank r = 3. Then, the AWGN b 
power is a 2 = 1. Finally, the theoretical covariance matrix of 
the total noise can be written as R = ^.^ UAU g + cr 2 I m 
with A = diag([6, 2,1]) and where JNR is the jamming to 
noise ratio, ypfipj is fixed at lOdB except for Fig. 4. 

In order to validate the spiked model as covariance matrix 
model, we visualize a zoom of the experimental PDF of 
the eigenvalues of our data without target in Fig. 3 over 
5x 10 4 Monte-Carlo iterations. We observe a Marcenko-Pastur 
law around 1 (eigenvalues of the white noise) and Gaussian 
distributions for the eigenvalues of the jamming, which is 
consistent to the CLT for the spiked model proved in [22]. The 
spiked model is consequently relevant for our data model. 

Moreover, in order to verified that the spiked model 



Eigenvalues 


where Eq.(32), Proposition 1 and the continuous mapping 
theorem were used to prove Eq.(41). One can observe that, 
although the traditional estimator of plr is consistent when 
K oo with a fixed to, it is no more consistent when 
to, K —> oo at the same rate. It is also important to underline 
that the new convergence result leads to a more precise 
approximation of /5 lr than previous works [16]. Indeed, [16] 
proposes an approximation dependent on K and the approxi¬ 
mation proposed here depends on K (and of course on c) as 
well as on the parameter 0. 

V. Simulations 

A. Parameters 

As an illustration of the interest of the application of RMT 
in filtering, the jamming application is chosen. The purpose of 
this application is to detect a target thanks to a ULA composed 
of to sensors despite the presence of jamming. The response 


Fig. 3. Zoom of the experimental PDF of jamming plus noise data with 

c = 0.2 and MM = lOdB. 
tr(A) 

is realistic in terms of separation condition. Fig. 4 shows 
(u r — y/c) as a function of in dB. This figure will be 
the same for all to and K at a fixed ratio. We recall that, 
in order to satisfy the separation condition, one should have 
u> r — i Jc > 0. Consequently, we gladly observe that it is 
satisfied for jyyxy > 4dB for the majority of c even c > 2. 
Indeed, in practice, if the jyyxy is lower, the jamming will not 
have any effects on the performance. 

B. Performance of filters 

We now observe the performances of filters through the 
SINR loss. We are first interested in the validation of the 
convergence of /Slr in Eq.(41) as to, A —► oo at the same 
rate. This convergence is validated and presented in Fig. 5 in 
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Fig. 4. Separation condition {uj r — y/c) of the spiked model for the lowest 
non-unit eigenvalue as a function of the ratio in dB. 


terms of MSE over 10 3 realizations with c = 3 for an AoA of 
the target (9 = 50°) and an AoA of the jamming (9 = 20°). 



Fig. 6 shows the visualization of Eq.(14) (blue line 
with stars), Eq.(15) (blue dashed line), the right side of 
the convergence in Eq.(41) (green line with circles) and the 
approximation E[/5 lr] — 1 — Jz introduced by [16] (black 
line) as a function of K when the target is near from the 
jamming, i.e. 9 = 20.5°. We observe that the spiked model 
and the RMT helps us to obtain a better estimation of ELolrI 
than the estimation E[/5 lr] — 1 — as the curve of 
has the same behavior as the curve of p lr. Then, similarly, 
the same equations are visualized as a function of 9 in Fig. 7 
with K = 2 r. We observe that, unlike the estimation 1 — r/K, 
the RMT with the spiked model permits us to obtain a better 
estimation of E[/5 lr] as a function of 9 and consequently a 
better approximation of its behavior. Thus, it permits to predict 
the parameter 9 value corresponding to the performance break 
(here around 21.1°). 



Fig. 6. Visualization of Eq.(14) (red line with squares), Eq.(15) (blue dashed 
line), the right side of the convergence in Eq.(41) (green line with circles) 
and the traditional estimation of E[plr] (black line) as a function of K (over 
10 3 realizations) with = lOdB, m = 100 and 6 = 20.5°. 



Fig. 7. Visualization of Eq.(14) (red line with squares), Eq.(15) (blue dashed 
line), the right side of the convergence in Eq.(41) (green line with circles) 
and the traditional estimation of E[/>| | ; ] (black line) as a function of 6 (over 
10 3 realizations) with z^ttI = lOdB, m = 100 and K = 2r. 


VI. Conclusion 

In this paper, we proposed new results in random matrix 
theory with a specific covariance matrix model fitted to our 
data model: the spiked model. Based on this, we studied the 
convergence of the traditional estimators of the SINR loss 
in their full rank and low rank version when the number of 
secondary data I\ —> oo with a fixed data dimension m and 
when to, K —>• oo at the same rate c = m/K. We observed 
that the full rank version is consistent in the two regimes. 
However, the low rank version is consistent when K —> oo 
with a fixed to but is not consistent when m, K —> oo at the 
same rate c. Finally, we applied these results to a jamming 
application. We first observed that the experimental probability 
density function of the eigenvalue of the covariance matrix of 
jamming data is relevant with the probability density function 
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of the spiked model. Then, we validated the convergence of 
the SINR loss in its low rank version and we observed that 
random matrix theory and more precisely the spiked model 
better evaluate the asymptotic performances of the low rank 
SINR loss corresponding to the adaptive LR filter, especially 
when the steering vector parameter is close to the jamming 
one and contrary to previous works. Moreover, it permits to 
predict the steering vector parameter value corresponding to 
the performance break. 

VII. Appendix 

The proof is decomposed as follows. We first develop the 
structured QF as a sum of simple QFs and base structured 
QF (Subsec. VII-A). In a second time, we formulate the 
base structured QF as a complex integral (Subsec. VII-B) 
and split it into several integrals (Subsec. VII-C). Then, we 
determine the deterministic complex integral equivalent of the 
base structured QF (Subsec. VII-D) and its formal expression 
(Subsec. VII-E). Finally, we use this result to determine the 
convergence of the structured QF in the large dimensional 
regime (Subsec. VII-F). The regime of convergences in the 
Appendix, if not precised, is m, K —» oo at a fixed ratio c. 


A. Development of the structured QF 

Let s i and s 2 be two deterministic complex vectors and 
B be a m x to deterministic complex matrix with uniformly 
bounded spectral norm for all to. In order to obtain the 
convergence of the structured QF s^II^BII^s^ one can 
rewrite, using the notations of Eq.(23) and the spiked model, 
ri^ = rir +1 = U r+ iUf +1 = I m - £l=i rii where 
U r+ i = [u r+1 ,--- ,u m ], ITi = UiU f, Vi G [1,r] and Ui 
are the eigenvectors of the SCM. We recall that r is fixed for 
all to, i.e. r/m —> 0 + . Thus, one can develop the structured 
QF as : 


H tt±-dtt-L H 

Sj Il c BII C S2=S X 


I m - B I m — ft, )s 2 (42) 


= Si Bs2 — Si ^2 riiBsa-sf B ^ n,S2 
2 = 1 2 = 1 

r r 

+ 8 ?'En i B'£n i 82 (43) 

2=1 2 = 1 

r 

= S fB S2 -]T (sfn 8 Bs 2 + sfBri,s 2 ) 


I 'y .s i TT I j fT , .s V .s ] nS Bn} 2 s 2 (44) 


ii=i 


j 1,32=1 

313^32 


B. Formulation of the base structured QF as a complex 
integral 

Remarking that Eq.(44) is a sum of simple QFs and base 
structured QFs, we first focus on the convergence of the base 
structured QF r)(j 1 ,j 2 ) = s^tl jl Bil h s 2 , {j 1 J 2 } G [l,r] 2 . 
Let us now formulate the base structured QF as a complex 
integral. 

Proposition 2: Let B be a to x to deterministic complex 
matrix with a uniformly bounded spectral norm for all to. 


Then, under (Asl-As5, As6.S) and the spiked model, Vji, j 2 G 
[l,r + 1], if t)(ji,j 2 ) = sfri^Bf^sa: 

31 32 

XB^R — Z 2 I m ) S 2 dz\dZ 2 (45) 


Proof: If j 1 fz j 2 , it can be easily shown that j 2 ) can 
be expressed as the following Cauchy integral: 


A = 


(2m 


sf(R — Zllm) ^(R — Z 2 Im) 1 s 2 dz\dz 2 


(46) 


where Cj in a negatively oriented contour encompassing the 

eigenvalues of R corresponding to the j -th eigenvalue of R 
and z 1 and z 2 are independent variables. Indeed, let G(zQ = 
(R-Zfclm)” 1 = (IXX^zi)- 1 with k G {1,2}. Thus: 


A = 


(2m 


(2m 


sfG(zi)B(R-z 2 I,n) l S 2 dz\dz 2 (47) 


sfG(zi)Bl y^X n u n u^-z 2 lr 


(2Z7r) 2 


xs 2 dz 2 dzi (48) 
f; a?G ^ l)B ^ 82 dz2dz 1 (49) 


A„ - Z2 


1 

2l7T 


Jc il n=l JC j 2 

where fi 2} (z 2 ) = °? G ^Bu n ufs 2 


(50) 


fr?\z 2 ), one observe that fif 1 (z 2 ) has a single simple pole 
\ n which is encompassed by C“ for the indexes n G A ij 2 
where A ij 2 is the set of indexes corresponding to the j 2 -th 
eigenvalue of R. Consequently, from complex analysis: 


X n —Z2 
r(2), 


From the expression 


A = 


1 

2*7T 


E/ s?G{ * l)B ' an ' a ” S2 dz 2 dz 1 


31 n£Mj 2 


2in T r - 


A „ — Z 2 


(51) 

(52) 


l/ c - E 2 f^d^dz, 

St rieM j2 JC j 2 

E [-Res(/^ 2) (z 2 ),A„)] dz! (53) 


31 neMj 2 


where Res (^fn\z 2 ), A„j is the residue of f$?\z 2 ) at X n 
Thus, using the residue theorem and residue calculus: 

1 


A = 


2in 


E 


3i n€-Mj 2 


- lim (z 2 - Xn)fn\z 2 ) 

z 2 — ^X n 


E _ lin { \ - 


~\7l n£A4j 2 


Z 2—*X. 


*2, 


dz x (54) 
s^G(zi)Bu n u^s 2 


An - Z 2 


dz\ 

(55) 


E <£ V lim (sfG(zi)Btt n ufs 2 } dzi (56) 

JC7 z 2 3\ n 1 > 


Ji n£Mj 2 
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2 * 7 T 


2m 


A= 2^T- S?G {Zl)B £ U n Un S2dZl 


n£A4 


sf (R - Zil„ 


32 

-1 


Bn,- 2 s 2 <fo 




^ ^ \ n U n U n -2-11???. I Bn i2 «2d*l(59) 


= 2 toi-Y, 


^1 'U'n'U'n Bn i 2 S 2 

An 2 l 




m r 


n=l 

f (i) 


A = 


£ ini- f - 1){Zl)dZl 

neM h JC j 1 

- ReS (/n^fcO.An) 

- 5 Z lin - («1 - ArO/^Ol) 


nEAij 


2 l —^ Atj, 


^ ^ 1 'U'n'U'n BU j2 S 2 

An 2 i 


yf lim < (An - zi 

_ . . 21 —>An I 

nSAIji n f 

yf lim < sfu„uf Bn j2 s 2 } 

. 21 —>\ n ^ ' 
n€A4j 1 1 


sf 


s(n.) = 


2 * 7 T 


c - n=l 


g(n« 

An 2 ^ 


-dz 


v{ji,h) = 


(2i 


sf 


(* 


2^1 In 


xB 221 m ) S2dZldZ2 


with 


(57) 

(58) 


Q(*fc) = YY h — « fe I m ) _1 


(72) 


(60) 

(61) 


where fn^Zi) = Sl u """_^ Ja82 . Similarly, fn\zi) has a 

single simple pole \ n which is encompassed by C“ for the 
indexes n £ M.j 1 ■ Thus: 


H(« fc ) = l m + « fe n(I m + fi) 1 U' ff Q( 2 : i: )U (73) 

Then, replacing (R — Zfcl m ) _1 by Eq.(71) in Eq.(45) and 
developing the obtained result, one obtains: 

VU 1 J 2 ) = 'f_ j> _ s?E(z 1 )BE(z2)s 2 dz 1 dz2 

C 01 C 32 

-( 2 wl- l- K<»)H(n)-C a (n)B] 

31 32 

xE(z 2 )s 2 dzidz 2 

-J2Wf c - i- [BC?Mfl(»)-‘ 

Jl J 2 

x §2(22)] dz\dz2 
ef( 2 i)H( 2 i)- 1 C 2 («i)B 


(2i7r) j c - j c - 


(62) 


xCf ( 22 )H( 2 2 ) 1 e 2 (z 2 )dzidz 2 

(74) 


D\ — D 2 — D 3 + D 4 

(75) 

(63) 

with 




E(z) 

= (I m + P)- 1/2 Q( a )(I m + P)- 1/2 

(76) 

(64) 

ef (z) 

= af(I m +P)- 1/2 zQ(z)U 

(77) 


C 2 (z) 

= n(i m + n)- 1 u if Q(^)(i m + P )- 1/2 

(78) 

(65) 

C ?{z) 

= (I m + P)- 1/2 ^Q(2)U 

(79) 

e 2 (z) 

= n(I m + f2) _1 U ir Q(^)(I m + P)“ 1/2 s 2 

(80) 


( 66 ) 


y U„Un BI[j 2 S 2 = sf nj 1 Bn j2 S 2 (67) 


Consequently, fj{j 1 ,j 2 ) = A for j- t j 2 . 

Then, if j-\ = ]■> = j and using the same arguments as 
previously, one has: 

h tx rjrr 1 I H ST u n u^Bu n u^ 

Si njBIIjS 2 = 7 — <p s 1 > -^- S 2 dz ( 68 ) 

Jc J n—1 A„ - S 

However, the remaining of the proof is based on the fact that 
the resolvent G(z) of the SCM can be found in the complex 
integral, which is not the case in the previous equation. 
Consequently, noticing that: 


D. Determination of the deterministic complex integral equiv¬ 
alent 

The convergence of the terms Di to I) \ has now to be 
studied. Some of them will tend to 0 and the remainder of the 
terms will tend to a deterministic integral equivalent. 

Proposition 3: Let B be a to x to deterministic complex 
matrix with a uniformly bounded spectral norm for all to. 
Then, under (Asl-As5, As 6 .S) and the spiked model, Vji. j 2 £ 
[l,r + 1], fjijuh) - VU 1 J 2 ) -^4 0 with 

v{jl,h) = (2^)2 j _ j _ ef (2l)H(2i) _ 1 C2(2l) 

xBCf (z2)B(z2)~ 1 e 2 (z2)dzidz2 ( 81 ) 


(69) 


where 7 ^ is a deterministic negatively oriented circle only 
enclosing 77 (cf. Eq.(26)) and 


where g(.) is a functional, Eq.( 68 ) is equivalent to Eq.(45). 
As a consequence, Vji, j 2 £ [1, r + 1]: 


H(*) 

— I m 3“ zb m (z)S2(I Tn + H) 

(82) 

ef (z) 

= Zb m (z)s?(lm + P)- 1/2 U 

(83) 

C 2 (*) 

= 6 m (*)n(i m + n)- 1 u ff (i m + p)- 1 / 2 

(84) 

cfw 

= zb m (z)(Im + P)~ 1/2 U 

(85) 

e 2 (z) 

= b m (z)S2( I m + n)- 1 ^ (I m + P)- 1 / 2 S 2 

(86) 


(70) 


C. Development of the complex integral 

Next, one want to split the previous line integral into 
several line integrals where some of them will tend to 0. Thus, 
from [ 22 ], with k £ { 1 , 2 }, one can write: 

(R-^Im)- 1 = (I m + P) -1/2 [Q(zfc) — ZkQ(zk)U 
xH(z fc ) _1 fi(Ir + n)- 1 U H Q(z l )] 
x(I m + P)“ 1/2 (71) 


Proof: We first recall that we are interested in the indexes 
ji , j 2 € [1, r]. Then, the function E(z) in D\, D 2 and D 3 
can be rewritten as: 


E(z) = (R - z(l m + P)) 


,H 


A ra — z( 1 + ui n ) 

A„ 


(87) 


Thus, E(^i) (resp. ~E(z 2 )) has a single simple pole .TO” ^ A„ 
when ui n 0, i.e. n £ [1, r] ((As5, As 6 .S) are verified and 
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f[x) —► f(x), with probability one for all large to, K at a 
fixed ratio c). As a consequence, V4,4 G [l,r], Cf (resp. 
Cf) does not encompass E(zi) (resp. E(z 2 ))- Thus, Di = 
D 2 = D 3 = 0 and: 

flUuh) = 7^-t^ / <f _ ef(2i)H(2i)~*C 2 (zi)B 

( tt) Jc H Jc j2 

xC?(z 2 )U(z 2 )- 1 e 2 (z 2 )dz 1 dz 2 (88) 

We will then determine a deterministic equivalent of 
Eq.( 88 ), i.e. its convergence in the large dimensional regime 
from lemma 5 of [32]: 

sup\\V H (Q(z) - & m (z)I ro )U|| ^4 0 (89) 

zeC m,K—too 

m/ K —^c<oo 


with 

g = ||r /_ ef (zi)H(zi)- 1 C 2 ( 2 i)d«i (103) 

A i 

in order to determine g in a first time. 

We recall that, in our case, wi > • • • > u r > ui r+ \ = 0. After 
an eigendecomposition of e( / (z\) and C 2 (zi) and, noticing 
from [ 22 ] that: 


H(z) _1 = diag 


l+zb m (z) 1+ l, 1 ’ ’ l+zb m (z) 


r +1 

^ 1 + Zb m ( 2 ) 3 ^ 


(104) 

(105) 


with 


with C a closed contour of C. Indeed, one can notice that: 

H(z) = lm + zSl(I m + fi)- 1 [U ff Q( 2 )U] (90) 

ef (z) = sf(I m + P)“ 1/2 zU[U h Q(z)U] (91) 

C 2 ( 2 ) = n(I m + fi)- 1 [U H Q(z)U]U"(I m + P)- 1/2 (92) 

Cf ( 2 ) = (I m + P)- 1/2 zU[U"Q(z)U] (93) 

e 2 (z) = n(I m + n)- 1 [U H Q(z)U]U"(I„ l + P)- 1/2 s 2 (94) 
Thus, from Eq.(89), one obtains: 

H(z) ^ H(z) = I m + zb m (z)fl(l Tn + n ) _1 (95) 

m,K—>o o 

m / K —tc<oo 

ef(z) 254 e f(z) = z5 m (z)sf(I m + P)- 1 / 2 U (96) 

m,K—t 00 
m / K —^c< 00 

C 2 (z) 24 . C 2 (z) = 5 m (z)fi(I m + n) _ 1 U H (I m + P )~ 1/2 (97) 

m,K —>-oo 
m/K—>-c< 00 

Cf(z) 24 Cf(z) = zb m (z)(l m + P)- 1 / 2 U (98) 

m,K—t 00 

m/K-)-c< 00 

e 2 (z) 24 e 2 (z) = 6 m (z)n(I m + n)- 1 U H (I m + P)- 1 / 2 s 2 (99) 
m,K—>oo 

m / K —>■ c< 00 

As a result, fj(ji,j 2 ) - t?(4,4) ^ 0 with 

i?( 4 , 72 ) = ( 2 ^ 2 j_ f _ ef ( 2 i)H(zi) _1 C 2 (zi) 

xBCf (z 2 )H(z 2 ) _1 e 2 (z 2 )dzidz 2 (100) 

where 7 ^ is a deterministic negatively oriented circle only 
enclosing Tj (cf. Eq.(26)). 

£. Determination of the expression of the deterministic equiv¬ 
alent 

Let us now find the expression of the deterministic equiva¬ 
lent 77 ( 4 , 72 ) as a function of the eigenvalues and eigenvectors 
of the covariance matrix R. 

Proposition 4: Let B be a m x to deterministic complex 
matrix with a uniformly bounded spectral norm for all to. 
Then, under (Asl-As5, As6.S) and the spiked model, 

77(71,72) = XilX02 S 1 n jl 1111)2 s 2 
with x i = 1 +1 1 - 1 and { 31 ,h} e [l,r] 2 . 

Proof: We first rewrite Eq.(100) as: 

77 ( 71 , 72 ) = ^ f - g X BC ?(^ 2 )H(z 2 ) _1 e 2 (z 2 )dz 2 

4 2 


( 101 ) 


( 102 ) 


Xi = 

one obtains: 


®Kl+...+«I- 1 


Ik, 


W + l 1 • • • 1 /N “r-(-l 


G C mXm (106) 


e l 


(zi)H(z 1 ) _1 C 2 (zi) = s? £ TTlffer 


(107) 


+"i 


Thus, 


r+1 _ 

H _ 2 l b m( 2 l) 


24 f- 1 ^ (7+“l ) 2 ^+zibm(z 1 )jff- 


^j—dzi (108) 


4, 1=1 


r+1 


444*44 


Z=1 

r+1 


2 1 ^rn ( 2 1) 

2*71 J - l + z ll>m(zi)x+ 


dz, (109) 


V 4—sfn,— cf 1+ J lb ™ (zi) —dzi (110) 

*+ 1+ “‘ 1 2m L- i±2i+zi6 m (*l) 


From [22], iif 1 + z-| h nl (z ,) = 0 only for Zi = 74 and 
*i 0+I z x =t h 4 0, 4 G [l,r]. Hence, 1 + 4 ', has a 

single simple pole at r 7l , 4 G [l,r], As a consequence, with 

7i(z) = zb m (z) 


1+ “7T Sl njl 2i7T / - +Z1 6 m(zi ) 


vT— sfn,-, 

i+^1 J- 


1 TT 


1 + 07 


1 + CJ 


1 H TT 

Sl Iljj 


—Res 


2i2A41)- dzi (H1) 

01 '"01 

h{zi)b m (zi) 


l + 7 / x ’ J1 

—^+/i(zi) 


(Tji) 

V / Ei=r i:L 


h’ijj 1 ) 


( 112 ) 


(113) 


(114) 


from residue calculus and where (.)' = 4+ Then, observing 


dz 

that , T 1 , = as +4 il _|_ hfrjf) = 0, one finally has 


1 _ t + hfrji) 1 +^jT 


„ _ „ 77 TT 1 + /l ( 4 l) M 4 l)M 4 l) 

8 - — n “~h(tP) h^T~ (115) 
, + (116) 

Then, similarly, we write Eq.(102) as: 

» 7 ( 7 'i, 72 ) = gB x g (117) 
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with 


g = 2^/_ Cf(, 2 )H( 22 )- 1 e 2 (, 2)d , 2 (118) 


g 


n,, s 2 


Thus, similarly as with g, one deduces that: 

(l + /t(r J2 ))fe m (r J2 ) 1 

h'(r j2 ) 

As a result: 

Vtiijz) = t(Tn)t(Tj 2 )s? TliiBn^ 

with 

(1 + h(Tj))b m (Tj ) 

h'{rj) 

Finally, the last step consists in expressing £ ( Tj ) as a 
function of or,. Using Corollary 2 from [22], one expresses 

£(Tj) as: 


£(Tj) = 


(119) 


( 120 ) 


( 121 ) 


1 — cu> 

£(tj) = Xj = —— 

1 + CLJ. 


-2 

J 

-1 


( 122 ) 


As a consequence. 


vUi,h) = XjiXj 2 s 1 Yl n BIlj 2 S2 (123) 

m,K —>oo 
m/K—*c< oo 


with {ir,j 2 } G II,r] 2 . 


E Convergence of the structured QF 

From the development of the structured QF, we recall that 
the convergences of the simple QFs sf n,B.s 2 and sf BIIjS 2 > 
V* £ [l,r] can be easily determined from [22]: 


sf fliBs 2 

a.s^ 

Xisfn,Bs 2 

(124) 


m,K—¥ oo 

m/K—)-c< oc 



sf BfUs 2 

a.s^ 

Xisf BIIiS 2 

(125) 


m,K—^ oo 
m/K —tc<oc 




where \i is defined as in Section III.C. 

Then, also using Eq.(123) in Eq.(44), one obtains: 

H |-r _L t-j r _L a.. s. ff -| j 

n c Bn c s 2 —>• s 1 Bs 2 

m,K—>oo 
m/K->c<oo 

r 

-y: (sf xin,Bs 2 + sf B Xi n iS2 ) 

i=1 
r 

+5Z n ' B x, n.;i«2 

Jl = l 
r 

+ Y1 S ' Ail n ji BX.,2 n ;, 2 s 2 (126) 

31^02 
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